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Abstract 

A generalized version of the nonequilibrium linear Glauber model with q states in d 
dimensions is introduced and analyzed. The model is fully symmetric, its dynamics be- 
ing invariant under all permutations of the q states. Exact expressions for the two-time 
autocorrelation and response functions on a d-dimensional lattice are obtained. In the 
stationary regime, the fluctuation-dissipation theorem holds, while in the transient the 
aging is observed with the fluctuation-dissipation ratio leading to the value predicted 
for the linear Glauber model. 

PACS numbers: 02.50.Ey, 05.50.+q, 05.70.Ln, 75.10.Hk 

1 Introduction 

The model introduced originally by Glauber [I] is denned on a one dimensional lattice 
where a one-spin-flip Markovian stochastic process takes place. It simulates the dynamics 
of a ferromagnetic Ising chain with first-neighbor interaction, and the stationary state is 
described by the Gibbs measure associated with the Ising Hamiltonian as a consequence 
of obeying detailed balance. In stationary regime at non-zero temperatures, the two-time 
functions like the autocorrelation and response function, are time-translationally invariant 
and connected through the fluctuation-dissipation theorem; furthermore, the model displays 
a transient where aging is observed [2J. At zero temperature, when the system becomes 
critical, the fluctuation-dissipation ratio (which is related to the effective temperature of the 
system[3j) assumes a nontrivial value Xoo = 1/2. 

When the model is extended to higher dimensions, and called Glauber model, no analyt- 
ical solution is available due to the nonlinear structure of the transition rate. Nevertheless, 
the linearized version of the Glauber model, proposed a few years ago [4], can be treated 
by analytical tools in any dimension. The linear Glauber model can be seen as the voter 
model with noise, which displays a disordered (paramagnetic) phase only; in the absence 
of noise, however, the system becomes critical. The dynamics of linear Glauber model was 
investigated by some of the authors in a previous paper [5]. This model is microscopically 
irreversible, that is, it does not obey detailed balance in the stationary state, for dimension 
d > 1. Moreover, it has both the stationary and aging regimes, with a nontrivial fluctuation- 
dissipation ratio Xqo = 1/2 in the later regime as in the usual one- dimensional Glauber 
case. 
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The fluctuation-dissipation relation was usually conceived for systems that obey detailed 
balance [61 [7], and has been generalized to include non-stationary regimes by the introduc- 
tion of an effective temperature measuring the violation of fluctuation-dissipation theorem 
(jS], and references therein). Many works have confirmed this phenomenon for several mod- 
els [9j QUI [HI H21 [13], HU H5]. Recent progress has suggested that it can also be invoked 
for nonequilibrium models [5j [16], which does not have an associated known Hamiltonian 
(see also this issue in the context of kinetically constrained models [T7J HE]). In [IB] , the 
fluctuation-dissipation relations were analyzed in a general class of models that exhibit up- 
down symmetry which does not obey detailed balance. 

The above mentioned works are closely related to the question of universality in out- 
of-equilibrium processes [19] . In equilibrium statistical mechanics, it is widely known that 
the critical behaviour of a system is governed by the fixed point of the renormalization 
transformation, and it turns out that only a few characteristics of the model are relevant to 
determine its universality class. 

In out-of-equilibrium dynamics, some of the problems related to universality may be 
addressed by the generalized version of the fluctuation-dissipation theorem 



where R(t,t') and C(t,t') are the response function and autocorrelation, respectively (see 
[2TJ] for some recent results). The typical experimental situation under consideration is a 
quench from a completely disordered state, which in reversible systems corresponds to a 
high temperature state, to the critical point. The usual fluctuation-dissipation relation is 
verified when the fluctuation-dissipation ratio X(t,t') equals the unity. It was conjectured 
that X(t, t') would depend functionally on C(t, t') only [2*1] , but renormalization group anal- 
ysis [22J and numerical calculations |23J indicated that the fluctuation-dissipation ratio is a 
function of t/t'. Furthermore, scaling arguments were casted to propose the asymptotic be- 
haviour for the autocorrelation and response function [10] at the critical temperature. This 
result suggested that the quantity 



is universal due to its dependence to dynamical exponents and the ratio of autocorrelation 
and response amplitudes, which are conjectured to be universal [2]. 

In this paper, we introduce and analyze a nonequilibrium lattice model, which is a gen- 
eralization of the nonequilibrium linear Glauber model to more than two states, which we 
call linear g-state model. The model can be understood as the linearized version of the 
dynamics associated to the equilibrium g-state Potts model. The dynamics of the model, 
as is the case of any dynamics of the Potts model, is invariant under the permutation of 
any two states. The model can also be understood as a g-state voter model with noise. In 
this interpretation, a group of individuals are called to vote in one of q candidates. A voter 
changes his opinion by choosing randomly a neighbor individual and adopting the neighbor's 
opinion with probability /i and remaining with his opinion with probability 1 — /i, the noise. 
Without noise, it reduces to the ordinary g-state voter model [21]. Similarly to the linear 
Glauber model, the present linear g-state model displays a paramagnetic phase whenever 
< [A < 1 and becomes critical at /i = 1. 

The analysis of such model has two main aims. First, it addresses the question raised in 
a previous result [5] about the dynamical phenomena of aging and violation of fluctuation- 
dissipation for a class of systems that does not obey detailed balance - recall that these 
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problems were usually studied through models that are described by a Hamiltonian. Finally, 
there is an additional interest in considering a model with a more general symmetry in order 
to verify its influence on the (possibly) universal quantity cited above, since in equilibrium 
statistical physics, symmetry plays a major role in the critical behaviour. 

The equilibrium Potts model with q states has been used to describe experimentally 
systems that display a number of identical states or structures at low temperatures |25j such 
as the adsorption of noble gases on graphite [25] . It has also been used to describe biological 
cell sorting [26J. The nonequilibrium model with many equivalent states such as the one 
studied here may be relevant in the description of systems where microscopic reversibility is 
not ensured like some biological phenomena [26 j. 

The layout of this paper is as follows. In section |2} the nonequilibrium linear g-state 
model is defined and many one-time functions are determined analytically The two-time 
functions are calculated in section |3l where the fluctuation-dissipation relations are carefully 
examined, and some dynamical exponents are calculated in section HI The summary of the 
main results and its discussions are found in the last section. 



2 The linear g-state model 

Consider a d- dimensional hypercubic lattice with N = L d sites and periodic boundary con- 
ditions. To each site % there is a spin variable Oi that takes the values 0,1,..., q — 1. 
The time evolution is governed by a one-site dynamics in which the state of a given site i 
changes from Oi to o[ = Oi + a modulo q, where a is one of the q states, and the states of 
the other sites remain unchanged. The possible transitions are then the ones in which the 
state o = (oi, a"2, . . . , Oii ■ ■ ■ i on) changes to the state o l ' a = (<7i, 02, . . . , o\, . . . , on) where 
a,- = Oi + a modulo q. The corresponding transition rate is denoted by wf(o) and, for the 
nonequilibrium linear g-state model, is defined by 



w i{°) = ~ + ^iJ2°( a i + a ^^+s)^ ( 3 ) 

where the summation is over the nearest neighbors and 5(x, y) is the Kronecker delta, which 
equals 1 if x = y and otherwise and the parameter fj, takes values in the interval < /i < 1. 
The time evolution of the probability P(o, t) of finding the system at state o at time t is 
governed by the master equation 

iP(°, *)=££ [<V^)PK<-\ t) - wf(o)P(o, t)} , (4) 

where the summation in a extends over the q states. 

The probability of a spin at site j be at state, say 1, is given by (6(oj, 1)). Throughout 
this paper, the notation 

(A(o)) = Y,A(o)P(o,t) (5) 

will denote the average over spin configurations of the state function A(o). The equation of 
motion for (S(oj, 1)) can be written from the master equation (jl]) as 

j t (o{o v 1)} = ?—J± - (Sfa, 1)) + ^ J2(5(o j+S: 1)). (6) 
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It is also possible to describe the time evolution of (8((Tj, l)5(a k , 1)), 



j t (5{a h l)5(a k , 1)> = [^-j 1)) + (5(a k , 1)>] - 2(S(a j , l)5(a k , 1)> + 

5 

which is closely related to the pair correlation. In order to recover the results obtained by 
jl] and [5] for the linear Glauber model, it is necessary to connect 5(<jj, 1) to an Ising spin 
variable Sj, that takes the values —1 or +1, through the relation 



s 



2<y(<7 i ,i)-i. (8) 



As a last remark on the model, its irreversible property for d > 2 will be discussed. 
Consider, for instance, the four states shown in figured] on a square lattice (d — 2). Suppose 
that the system follows the sequence of states A, B, C and D and returns to the initial state 
A. If the interval between two successive states is At, the probability of occurrence of the 
sequence A — > B — > C — > D — > A can be calculated through the transition rate Wi (which is 
wf =0 (a)) as 

P(A B C D A) = P(A\D)P{D\C)P(C\B)P(B\A)P{A) 

+^ ( l -^ + 2^\ . (9) 



\ q 2d ) \ q 2d) \ q 

This result is not necessarily equal to the probability of observing the reveresed sequence 
A^D^C^B^-A, which is 

P(A^D^C^B^A) = P(A\B)P(B\C)P(C\D)P(D\A)P(A) 



\ q 2d) \ q 2d) \ q 2d, 

The existence of a sequence of states that is not reversible implies that the system is irre- 
versible. A generalization of this result to higher dimensions is obtained by, for instance, 
filling the sites created by the introduction of more dimensions with spins a where a ^ 0, 1 
(this is the situation when q > 2; the case q = 2 was already discussed in [5]). The argument 
above can also be invoked to show that the model is reversible in the one-dimensional case. 



2.1 Site magnetization 

The definition of the site magnetization will be guided by some constraints. At a fully 
ordered state, where the spin is at state, say <x,- = 1, one should have (5(<Tj, 1)) = 1, while 
at disordered state, where the spin is at any one of the q state with equal probability, the 
condition (<5(oj, 1)) = 1/q should be satisfied. This leads to a natural definition of an ordered 
parameter Y^j rrij/N, where the site magnetization is defined by 

m j (t) = -±-(5{* j ,l))-^ 1 . (11) 
q — 1 q — 1 
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Figure 1: A possible irreversible sequence for the bidimensional (/-state model (q > 2). 

The above definition leads to 

■gmM = -m.it) + |j (12) 

s 

which is the time evolution of the site magnetization. 

The equation ffl2l can be solved by the usual methods by introducing, for instance, the 
Fourier transform 

<(t) = E^(*) e ^ ( 13 ) 

i 

and its inverse 

^)4E<(*)e*. ( 14 ) 



N v 



where the summation in p is over the sites of the the first Brillouin zone, in which each 
component of the vector p takes values inside the interval between — it and 7T. The solution 
of the differential equation ffl2l is then 

m j (t) = Y / r j -e(t-t')m e (t'), (15) 

l 



where 



r ^) = ^E^ W " (p)i , (16) 
i v p 

and 

/» = l-4Ecosp,, (17) 
a i=i 

for a d- dimensional hypercubic lattice with coordination number z = 2d. For a homogeneous 
initial condition rrij(t') = m° for any j, it is straightforward that the site magnetization 
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is constant {rrij(t) = m°) for /i = 1. On the other hand, the condition /i ^ 1 implies 
rrij(t) = m°e~( 1 ~^( t ~ t \ which means that the magnetization decays to zero for sufficiently 
long time. The time correlation length r, defined by m = m°e~ t l T ', diverges as r ~ (1 — , 
from which vn = 1. Moreover, the static magnetization M = \im t ^ oc J2j m j(t)/N is always 
zero for // ^ 1 and is a non-zero constant (if m° ^ 0) at criticality (/i = 1); this jump in the 
magnetization implies — 0. 



2.2 Pair correlation 



The definition of pair correlation qj t k{t) for spins at j and k will be gauged to obey the 
following requirements: (i) qjj{t) = 1 for any j, (ii) qj t k{t) = 0, j ^ k, for the paramagnetic 
state (iii) qj,k{t) — 1, j ^ k, for the ordered state. These conditions leads to a natural 
definition for the pair correlation, which is 



q 



- ^[mj(t) + m k (t) 



1 



q 



;is) 



and is the only one where the pair correlation qj,k{t) is linear to (5(<7j, l)5(<Jk, 1))- Its time 
evolution, 



dt 



8 

(g-2) [m,(t)+m fc (t) 



(19) 



is obtained from the master equation (j4])and from equation (J7J). 

From now on, it will be assumed that the pair correlation qj^if) depends on sites j and 
k through their difference r = j — k only (note that q r {t) = q- r {t))- Moreover, the system 
will be assumed to be in a random initial state (see [27], [28] and [29J for other possibilities), 
such that the evolution process can be understood as a quench from fi = l/qtoa.fi^l/q, 
which implies rrij(t) = for any t. This condition leads to 



d 
dt 



Qr(t) 



d 



(20) 



The above equation, which is valid for r / only, should be modified to comprise the case 
r = 0, for which qo{t) = 1. By using a previously introduced method [3], we write the 
equation fTSOj) in the form 



d 
dt 



q r {t) = -2q r (t) + ^J2<lr+6{t) + b{t)5 rfi , 

s 



(21) 



which is now valid for any r, including the case r = 0, provided b(t) be chosen to ensure 
qo(t) = 1 (note that if r ^ 0, then equation (12T|) recovers (120|) ). Formally, this means that 
the function b(t) should satisfy 



b(t) 



d 
~dt 



q (t) + 2q (t) - ii £ q s (t) = 2 - £ £ g , (t) 



(22) 
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Since the system was assumed to be in a completed disordered initial condition, then 
g r (0) = 5 r fi and the equation of motion for the pair correlation ff2T]) can be written as 



r,0) 



where 



dt e- st q r (t) 



is the Laplace transform of q r (t) (similar formula connects b L (s) and b(t)). 
The equation ( 123|) can be solved by introducing the Green function 



irp 



where / is defined in ( ITT)) , that satisfies 



iV^ s + 2/» : 



sGf (a, n) + 2G L r {s, A*) - ^ E G £fA A*) = *r,o- 

Hence, the solution of the non-homogeneous differential equation (1251) is computed as 
g r L ( S ) = £ G£_ r ,(s, fi) \l + b L (s)]5 r , fi = \l + b L (s)]Gi:(s } fi). 



(23) 



(24) 



(25) 



(26) 



(27) 



The function b L (s) is fixed remembering that the condition q (t) = 1, or qo(s) = 1/s, should 
be satisfied. It is easy to see that 



b L (s) 



- h 



(28) 



which implies 



sG%(s,ii)' 



(29) 



The stationary value for the pair correlation is obtained through the Laplace final value 
theorem 



q r (oo) = lim q r (t) = limsg r (s) 

t— >oo 



G r L (0,/i) 



In the same fashion, the stationary value for b{t) can be calculated as 

1 



(30) 



b{oo) = lim b{t) = \vcsisb (s) 



t— >oo 



Gf(0,A*)" 



(31) 
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2.3 Susceptibility 

In this work, the susceptibility is defined through the (spatial) variance 
X(t) = £ [<%o, 1)> - <%o, l))(5(a r , 1)> 

r 

(q-l 



Q 



X)?r(f), 



(32) 



where the random initial condition is being assumed. Starting from a disordered state 
{rrij(0) = for any j) and assuming \i ^ 1, the stationary susceptibility, 

q-l 1 



X(oo) = lim x(t) - N^L/n v 

is obtained by invoking the previous result (|30l) . In a hypercubic lattice, one has 



(33) 



G£(0,//)~ 



H^J r H) (1 -" )! 

1 



27TyU 



In (1 - ^) 



< d < 2, 



d = 2, 



(34) 



as /i ~ 1, and lim^i Gq (0,yu) < oo for <i > 2. Therefore, the stationary susceptibility is 



' q - l\ /2vr 



r i-| 



7T/i 



(1-m) -1 



9 J 

'q-A i 

v 2g2 J G L (0,/i) 



In(l-M)] 



< d < 2, 



d = 2, 



d > 2, 



(35) 



from which the exponent 7 is obtained: the susceptibility diverges algebraically with expo- 
nent d/2 for < d < 2 and 1 for d > 2 with logarithmic corrections for d = 2. 



3 Two-time autocorrelation and response functions 
3.1 Two-time functions 

The analytical form for the autocorrelation and response function will be determined in this 
subsection in order to analyze the stationary and aging dynamical regimes. The two-time 
autocorrelation is defined as 

C(t,t>) = lim 1 £ [<%,•(*), !)%#),!)> - (<%(*), (36) 

iV— »oo iV ■ L J 

with the two-time correlation 

(%^t),i)% i (tO ) i)) = EE%i(*) ) 1 ) p (^*i^*0%K* / ) ) i)^(^* / )> (37) 

<7 cr' 
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where P(cr, t\a', t') is the conditional probability of finding the configuration a at time t given 
the configuration o' at an earlier time if . Noting that 



(5(aj(t), 1)> = £ 6(a 3 (t), l)P(a, t\a', t') 



(3f 



with the condition at time t' being (8(<jj(t'),l)) = 5(a'j(t'),l), and invoking the definition 
(TTTT) and the equation ffl5l) . it is possible to show that 



C(t,t') 



lim 

7V-s>oo 



g-1 



(39) 



for the disordered initial condition. 

On the other hand, if one assumes an arbitrary initial condition, one has 



C(t,t') 



lim 



'^1 E r,(t - f)*(0 + (g " 1) J g " 2) ^i:^) 



g-1 



— ^m^m^t') 



(40) 



The evaluation of response function requires the presence of a (small) perturbation on the 
system. In the analysis of stochastic models, the introduction of an external field modifies 
the one-spin-flip rate to 



1-M 

. 4 , 



i) + ^ E i) + o{h) 



(41) 



where a Taylor's expansion was performed in the last step. 

Performing similar calculations of subsection 12.11 it is possible to show that 

= _ TO .( t) + Rj2m j+S (t) + y q b(t)h,(t), 



(42) 



assuming again disordered initial condition, when rrij(t) = O(hj). The solution of this 
differential equation, which can be obtained following the same previous ideas, is 



= ^E f dt'T^ k (t-t')h k (t')b(t'). 



(43) 



The above result, 
R(t,t') = 



is sufficient to evaluate the two-time response function 



1 5(8(a,(t),l)) 

N^ooN 1 - 8hj(t') 



g-1 

- Q , 



8m j{t) 



lim — 



HO 



g-1 



ro(t-t')KO- 



(44) 



It is worth to stress that the formula ( 1441) for the autoresponse function is obtained even 
assuming an arbitrary initial condition. 
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3.2 Stationary regime 



The stationary regime can be realized when both the waiting time (t') and observational 
time (t) grow with the constraint that r = t — t' > is fixed. In this limit, and assuming 
disordered initial condition, the autocorrelation, 

1 r d d n 

CM=c{T)= m^J^vM- (45) 

and the response function, 

flftO = RiT) = (^i) ^|L, (46, 

are functions of the time difference r only, and they are related to the usual form of the 
fluctuation-dissipation relation R(r) = dfC(r), as expected in a stationary regime. 



3.3 Aging regime 

The aging scenario can be seen when both the observational time (t) and waiting time (t') are 
made large without the difference r = t — t' being fixed. In the stationary regime, where the 
difference t — t — t' was fixed, the limit t' — > oo made the function q r (t') in ( 1391) and b(t') in 
(jUj) time- independent. This is not the case in the aging regime, where both autocorrelation 
and response function depend on t and t' independently. More precisely, the transient is 
observed if t 3> t', and this condition can be realized if the limit t — > oo is taken before the 
limit t' — > oo. Assuming disordered initial condition, and from previous results, it can be 
shown at criticality /i = 1 that in this regime the autocorrelation function behaves as 



C(t,t') 



q 2 



2I+ 1 sin f ^ 



7T 



~ < 



V q 2 J (t-t')lnt' 
' q - 1\ / d V 1 



< d < 2 



d = 2 



V g 2 ; \2txJ G l (0,1) 
and the response function is asymptotically equal to 



'q-l\ dSin(^) a d 

7T 



(t - 2 t' , d > 2 



R(t,t') 



1 



v 2q 2 

V 2q 2 ) (t-t')lnt' 
'q-l\ { d\* 1 



2q 2 



2nJ G L (0,r 



(t - 



< d < 2 



d = 2 



d > 2 



(47) 



(48) 
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The above results agree with the scaling C(t, t ; ) ~ t' f c (t/t' ) and R(t, t') ~ t' a f R (t/f) 
where a = b = (d — 2 + rj) / z (see table I), and fc and fn are scaling functions that behave 

as fc/ R (t/t') ~ (t/t') _A/z for t/t' ~ oo. 

In the aging regime, the fluctuation-dissipation theorem is not expected to hold anymore. 
The fluctuation-dissipation ratio 



X(M') 



R(t,t') 
d t ,C(t,t') 



(49) 



which measures the distance of the model to the stationary state (when X(t,t') = 1), has 
the following limit: 



X(oo,t') = lim X(t,t') 

t— too 



6(t')/2 



b{t>) - (1 - n) X (t') 

where b(t) and x(t) are given, respectively, by (1311 and (1551) . This result implies 

1, ^1, 



(50) 



lim 



1 

2' 



(51) 



which is identical to the Ising case [5]. 

The previous result has considered disordered initial condition, which assumes rrijit = 
0) = for every site j. If one starts from an arbitrary initial condition, it is possible to show 
that (now using equations fj40|) and (JHJ) 



X(ooA') 



lim 



6(0/2 



6(f) 



(52) 



which shows that a similar formula for X(oo,t') is obtained even for an arbitrary initial 
condition. The fluctuation-dissipation ratio is identical to 051)) ; other non-trivial values 
for this ratio (for non-zero magnetization as initial condition) can be found, for instance, in 
j, [28] and [29]. 



4 Dynamical exponents 

4.1 Dynamical exponent 9 

From the solution of equation (|T2|) . 

rrijit) =m°e- (1 -^, (53) 

one sees that at the critical point \i = 1 the magnetization is constant and does not vary 
with time. This implies the exponent 6, defined through mj(t) ~ m°t e [31] in the short- 
time regime, to be zero. It is possible also to calculate this exponent by means of the time 
correlation of the total magnetization [32J . 



11 



4.2 Dynamical exponent X/z 

At the critical point fi = 1, one may calculate the dynamical exponents A and z, defined 
through C(t, 0) ~ t~~ . From ( l39~I) . it is immediate that 



The q r (0) can be evaluated by invoking (|29|) and the Laplace initial value theorem 

g r (0) = lim + g r (t') = Hm sq^(s) = 5 T:0 . 



In the thermodynamic limit, this result implies 



C(t,0) 



q- 



(27T///d) 



(54) 



(55) 



(56) 



where Jq(x) is the modified Bessel function of order and the last passage is obtained in the 
asymptotic limit f > 1. The autocorrelation decays exponentially for /i ^ 1; nevertheless, if 
\i = 1, one sees that 

- = i (57) 



4.3 Dynamical exponent ( 

Another dynamical exponent of interest is (, defined through x(t) ~ t^. From the results 
obtained in subsection 12.21 the Laplace transform of the susceptibility can be written as 



X L (s) 



(5? 



The asymptotic behaviour of the (dynamical) susceptibility x{t) f° r large times corre- 
sponds to the Laplace anti-transform of x L { s ) when s ~ 0. In this regime, one can evaluate 
the Green function at criticality as 



Go(M)~< 



4tt) r ( 1 



27T 



-S3. 

2j 



which yields 



x(t) 



2 (q- 1 \ ^47rV 



1 



d V <? 2 J yd J r (I) r (l 

q-l\ t 



2tt 



g 2 / lnt' 



1 



W 

, q 2 M(M)' 



< d < 2, 

d = 2, 
d > 2, 

*5, 0<d<2, 

d = 2, 
d > 2, 



(59) 



(60) 
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showing that 



d 
2' 



< d < 2, 
d > 2, 



(61) 



with logarithmic corrections for d = 2. 
4.4 Dynamical exponent z 

The exponent z, defined by the behaviour of the correlation length £ ~ (1 — /i) z will be 
estimated through the spatial correlation, which can be casted as 



q r (t = oo) 



d 

d \ 5 



/i\ V 2 4 (1 — jj) 



I 12(1- n)d 
-z—K d-2 r\ 



G#(0,/i) V 2 W ^ dJ G o (0,//)r 
where K v (z) is the Macdonald's function, which behaves as 

' 2 l/ ~ 1 r(z/) 
— , \z\ < 1 and v ^ 0, 

~ <v In f-^ , |z| < 1 and v 0. 



, |z| > 1. 

Therefore, for large distances, the correlation decays exponentially as e~ r ^, where 



(62) 



(63) 



(64) 



is associated to the correlation length. On the other hand, when the system is near criticality 
in the sense that r <C £, one has 



g r (t = oo) 



r (^) (r/20 
In (r/2£) 



d-2 ' 



ln(v^J/2£)' 

d r(g) i 

» 4vrf/iG'o(0,/i)r 



d-2 ' 



< d < 2, 



d = 2, 



d > 2, 



(65) 



from which it is possible to see that the exponent r/[33j is equal to zero for d > 2 and is 
equal to 2 — d if < d < 2. From fl64"j) . the exponent z/j_ (defined by £ ~ (1 — /i) ^) is 
equal to 1/2, while the dynamical exponent z = is then 2 = 2 (since Un = 1, as seen 

in subsection 12. ip . 
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Table 1: Critical exponents for the nonequilibrium linear g-state model 



Exponent 


< d < 2 


d > 2 


c 


d/2 


1 


z 


2 


2 


A 


d 


d 


e 









Exponent < d < 2 d> 2 
P 0~ 

^1 1 1 

z/ ± 1/2 ~l/2 

7 ~ d/2 T~ 
77 2- d 



Collecting the previous results (see 12.11 12.31 14.21 14.31 and I4.4R . the table [U is obtained. 
One should recall that for d > 2, one has £ = 1 and 7 = 1 (the other exponents remain 
unchanged). These values for the exponents satisfy the relations (z = (d — 2/3 /v) and 
9z = (d — A). It is worth mentioning that these exponents are in agreement with the 
universality class of the voter model [3U [35] . 



5 Conclusions 

This paper has established a number of exact calculations for the dynamical and static 
behaviour of the d-dimensional nonequilibrium linear g-state lattice model, which is a gen- 
eralization of the nonequilibrium linear Glauber model. This model is fully symmetric in 
the sense that it is invariant under the permutation among all the q states, having the same 
symmetry of the equilibrium Potts model. Although the analytical form of many functions 
are now distinct from the linear Glauber model, many similarities were reported in this pa- 
per. The stationary and aging regimes were both characterized, with the usual fluctuation- 
dissipation relation satisfied in the former regime and violated in the later regime at the 
criticality /i = 1. The fluctuation-dissipation ratio indicates that the dynamical be- 
haviour of the present model is similar to the linear Glauber model with Ising spins, which 
is just a particular case (q = 2), and thus independent of the number of states q. 

When ji — 1, we recover the voter model with g-state and in this case the system 
finds itself in the critical state. From the results for the correlation functions it is possible 
to generalize a result already known for the case q = 2 about the stationary states. In 
one and two dimensions the only possible stationary states are the ones in which all the 
sites of the lattice are in one of the q absorbing states, which one depends on the initial 
configuration. This statements seems to be in contradiction with the result (153]) which says 
that the magnetization remains constant. To understand this it suffices to remember the 
meaning of choosing an initial condition with magnetization m . This means that one should 
consider several initial configurations whose average gives the magnetization m . Each one 
of these configuration will reach one of the q absorbing states. The averages over these 
absorbing states will give an average m which according to ( 153"]) should equal mo- In three 
or more dimensions there are other states stationary states besides the q absorbing states. 
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7 Appendix 

In this appendix (and also in the main text), the Landau notation was adopted: 

i) if f(x) = 0(g(x)) (assuming g(x) > 0), then there exists xq such that \f(x)\ < Ag(x) for 
some constant A if x > xq. 

ii) if f(x) = o(g(x)) (assuming g(x) > 0), then hm x -^ 0O f(x)/g(x) = 0. 
7.1 Dynamical susceptibility (case d = 2) 

The asymptotic behaviour of the dynamical susceptibility is calculated through the Laplace 
anti-transform 

1 rc+ioo 

Xit) = — dse st X L (s), (66) 

where c is real and larger than the real part of any pole of x L ( s )> given by (!58|) . Since the 
< d < 2 and d > 2 cases are simpler, the evaluation of dynamical susceptibility will be 
presented for d — 2 only, which implies 

^i^lrfe' (67) 

One should first consider the integral 



"j rc+ioo p s ^ 

Bit) = / ds — -. — - — -, (68) 

w 2-KiJc-ioo s -Ins ' K ' 



Ins)' 

where c and c are real and larger than the real part of any pole of the integrand. The 
function B{t) relates to x(t) through 



X (t) = 2vr ( q —\ B(t). (69) 



By using the contour shown in figure |2] and invoking the residue theorem, the integral 
( 168]) can be casted in the form 



B(t) = / dr^- 

Jo r h n 2 r + 



7T 2 



B a (t) + B h (t) + B e (t), (70) 



where 



tint , e 



B a{t) = -v, (71) 



(in 2 r + 

15 



7T 2 




Figure 2: Integration contour used in the evaluation of the function B(t) given in 



and 



B b (t) 



B c (t) 



-rt 



dr- 



^ln 2 r + 



7T 



OO g 

dr 



-rt 



r(ln 2 r + vr 2 

These three functions, B a (t), Bb(t) and B c (t) will be evaluated separately. 
Function B a (t) 



Since 



and 



then 



B a (t) 



i 

t hit 



dr- 



L 



i 

t hit 



dr- 



11 r (in 2 r + 7r 2 ^ 
-rt)' 



1+E 



m=l 



ml 



(in 2 r + 



E 



7T m=l 



m! 



1 



o 



r In r + 7r 2 m=1 



£ (rt)' 



< 



rln 2 r Vilnt 

1 fthxt i 

ar- 



ln £ Jo r In r 
i i 





lntln(t hit) 



B a {t) 



i 

i In i 



dr- 

o r 



(in 2 r + 

-Aln(t]nt) dy 



•{ arctan 

7T 



7T (y2 + 1) 

1 In (tint) 



o 



1 

L77 



7T 
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where the change of variable r — > e ny was performed in the second line. For t > 1, one has 



B a (t) 



1 



In (t In t) 



O 



1 



\n z tJ hit 



J_ + (\n{\nty 



ln 2 t 



Function B b {t) 

The function Bf,(t) has the following upper bound: 

1 



\B b (t)\ 



-it 



dr 



_L_ r ^ n 2 r _|_ 



7T 



< 



— , e rt 
ar 



lnt P ~ u 



du- 



ll 



^ e"Tni r^t du 2e In (hit) ^ In (\nt)" 



In t 



\n z t 



\n z t 



Function B c {t) 

The function B c (t) has the following upper bound: 

\B c (t)\ = 



oo g- 

dr 



-rt 



< 



1 f 00 . e" r * 



1 1 f 00 
dr- — < ^r-r-r / dre~ rt 



IT I lILL V 7T l int 

t t t 



1 



O 



1 



vr 2 tlnt Vtlnt, 
From ([TOD, (EZD, (J7HD and ([79]), one fin ally has 



which can be inserted in ( 1691) to yield 

' q -r 



X {t) ~ 2tt 



2tt 



9 



, dtfl(t)=27r(^— ^ 



9' 



eft 



1 



lnt In 2 * 



dt 



ln 2 t 



g 2 / lnt 



(77) 



(7c 



(79) 



(80) 



7.2 Spatial correlation function 

In the thermodynamic limit (N — > oo), the spatial correlation function can be casted as 



q r (t — > oo) 



GV(0,/i) 



d d p 



G (0, /x) G (0, fi) J[-*,*)* (2vr) d e + 2fi (l - J £? =1 cosp, 
as seen in subsection 12. 2[ and e = 2 (1 — //) . It is straightforward to see also that 

1 f d d p < x 



(82) 



q r (t — > oo) 



G (0,/i) ./[-7r,7r)<* (2tt) 



1 



n 



G (p,n)Jo e i=1 
17 



cos p, 



J3) 



where I v (z) is the modified Bessel function, which behaves asymptotically as (|_36j) 

l v {z ~ oo) ~ e -j=- (84) 

V Z7TZ 

Therefore, when the system approaches the critical point (e ~ + ), one has 



d 



d 

2 jo /*00 r 2 £d 



q r (t ^ oo) = -— — - - — e2 / dye 



-y- 







( d ^ 1 PC 1 , ^ 

= o — 7w7^ — V ^ — A ^l r /?)' I 85 ) 

V 27r M/ G o(0,/i) (r/g)~ s ~ 

where ^(-z) is the Macdonald's function (or modified Bessel function of the third kind) 

1 / Z\ U f°° z 2 7T 

Kv{z) = 2\2) 1 d V e ~ V ~ lirV ~ 1 > l«S*l<4. ( 86 ) 



and 



f=yg(l-rt-*- (87) 
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